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(— I ' Abstract. We discuss the possibility to identify anisotropic and/or inhomogeneous 

<*S |. cosmological models using type la supernova data. A search for correlations in current 

q | type la peak magnitudes over a large range of angular scales yields a null result. 

i—{ • However, the same analysis limited to supernovae at low redshift, shows a feeble 

anticorrelation at the 2a level at angular scales 9 « 40°. Upcoming data from, e.g., the 
SNLS (Supernova Legacy Survey) and the SDSS-II (SDSS: Sloan Digital Sky Survey) 



supernova searches will improve our limits on the size of - or possibly detect - possible 
correlations also at high redshift at the per cent level in the near future. With data 



(N 
> 

\Q ' from the proposed SNAP (SuperNova Acceleration Probe) satellite, we will be able to 

0^ ■ detect the induced correlations from gravitational lensing on type la peak magnitudes 

on scales less than a degree. 
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1. Introduction 



Apart from the identification of dark matter, the most dominant question in cosmology 
today is that of what is responsible for the apparent acceleration of the universal 
expansion. Assuming that the cause is a dominant energy component with negative 
pressure, current efforts are focused on establishing whether this dark energy (DE) 
component can be described using a cosmological constant (CC), i.e. an energy 
component with constant density and a fixed equation of state (EOS) given by p = up 
where w = —1, or whether dark energy is dynamical (DDE) with a varying EOS, 
w = w(z). 

Two main techniques are employed in these investigations: probing cosmological 
distances and probing the growth of cosmological structures, both of which are sensitive 
to the energy content of the universe. Distances are most robustly probed via type la 
supernovae (SNe la) [TJ [2j [3j 0] , the baryon acoustic oscillation (BAO) peak [SIEIE] and 
the last scattering surface of the cosmic microwave background (CMB) |8J. Structure 
growth is mainly probed by large galaxy surveys such as the 2dF [9] and the SDSS 
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[TU| [TT] but also via galaxy cluster counts [12] and weak lensing [T3]. Despite some 
claims to the contrary, there is a consensus that current data are perfectly consistent 
with having approximately 70% of the energy density in the universe in the form of a 

cc ng. 

However, cosmological distances only depend on the dark energy EOS in the form 
of a double integral and it is thus very difficult to investigate the time evolution of 
the EOS parameter, w(z). Though not demanding it, current data thus still allow for 
large deviations of w(z) from the CC value of w = — 1 |14j . Even if the limits on w(z) 
will improve in the future, detecting such an evolution remains notoriously difficult, 
especially if the deviation from the CC value turns out to be small. 

It is thus of interest to look for alternative ways to detect or constrain any behaviour 
of the DE that differs from that of a CC. One such possibility is to study not temporal 
evolution, but rather spatial variations in DE properties. The CC value w = — 1 is 
special in the sense that it not only gives a constant energy density in time, but also 
that it prevents any spatial clustering of DE. Any deviations from this value inevitably 
cause the DE to cluster, a clustering which - if detected - would refute the CC as the 
dominant energy component in the universe. 

The clustering properties of DDE in different models and scenarios have been 
investigated in several recent articles (see, e.g., [15] and references therein). For values 
close to w = —1, such clustering is expected to be weak and take place mainly 
on very large scales, larger than the current Hubble radius and thus inaccessible to 
observations. On smaller scales, it can be shown that the DE density will tend to 
be anticorrelated with the matter density [151 IB]- Although, again, the amount of 
clustering is expected to be small (on the order of 10~ 5 [1 + w]), it is nevertheless 
important to get observational confirmation that this is indeed the case, even if it is 
not expected from purely theoretical considerations. On a more speculative note, also 
in models where the apparent acceleration of the universe is explained in terms of the 
back-reaction of inhomogeneities in the universe, we expect to have spatial variations in, 
e.g., cosmological distances [17] . Furthermore, cosmological observations can be used 
to put constraints on the degree of centricity in spherically symmetric inhomogeneous 
models [HI [H] and on anisotropic DE models [HI EH E21 E3 [24] . 

One way to probe models with inhomogeneous DE or back-reaction is to look 
for anisotropies in the observed peak magnitudes of SNe la. In general, we expect 
cosmological distances to be correlated on scales similar to the clustering scale of the 
DE or the matter inhomogeneities responsible for the back- reaction effect. The detection 
of such correlations is demanding because of the intrinsic variation and observational 
uncertainties in SN la peak magnitudes. 

Since systematic effects connected to, e.g., the observational properties of different 
telescopes, calibration issues and details of the light-curve fitting procedure can induce 
correlations in SN la magnitudes, it is important to minimize these effects by using 
as homogeneous a data set as possible in the analysis. Also, physical effects such as 
peculiar motions [55] . gravitational lensing [26] and dust extinction [27] will introduce 
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correlations in the SN la data and need to be controlled and, if possible, corrected for. 

In Kolatt & Lahav [28] , an early data set consisting of a total of 79 SNe la combined 
from the Supernova Cosmology Project [2] and the High- 2 Supernova Search Team [1] 
was used to look for directional variations in the best fit cosmological parameters. The 
result is consistent with the expected statistical variations in a homogeneous universe. 
Gupta et al [29] used a more homogeneous set of 157 SNe la [30] to look for directional 
variations in the statistical scatter around the best fit cosmological model, again with a 
result consistent with null variations. A slightly different technique was used in Bochner 
[31] to look for anisotropic scatter in 172 SNe la [32] with similar results. Although 
current investigations have been mainly inconclusive, Schwartz & Weinhorst [33J find 
an offset in the best fit calibration of low z SNe la, between the north and south 
equatorial hemispheres at the 95% confidence level. Whether this hint of a north/south 
asymmetry is due to a statistical coincidence, observational systematic differences, or 
correlated peculiar motions or has a cosmological origin is not clear at the moment. 

In this paper we devise a general methodology for detecting angular correlations 
in SN la magnitude residuals. We apply this method to current SN la data as well as 
simulated upcoming data sets. 

In section [5J we describe the method for analysing SN la magnitude residuals for 
angular correlations and define the detection limit for such correlations. Two data sets 
from the literature are analysed for angular correlations in section [3j In section HI we 
present different SN la surveys and investigate the detection limits for each of them. In 
section [5l we assume a toy model correlation function and illustrate how future data 
can put limits on such a correlation. The paper is concluded in section [61 

2. Method 

In this paper we study angular correlations only, since the effect from possible DE 
inhomogeneities is larger in the transverse direction than in the radial direction. 
Moreover, a full three-dimensional study would require larger statistics than currently 
available. The data sets we investigate contain both nearby and distant SNe la, for 
which we have redshifts, z, observed peak magnitudes, m obs , with uncertainties, a m , as 
well as positions on the sky, right ascension, a, and declination, 5. 

2.1. The angular correlation function 

The magnitude residual of a SN la on the Hubble diagram is the difference between 
its observed magnitude and the magnitude predicted from its redshift in the best fit 
cosmology, 

5m = m obs - met . (1) 
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We define a normalized and scaled residual, X , by subtracting the mean and dividing 
by the uncertainty, 

x = Sm - \Hrn ^ 

Note that the o~ m include both observational errors and the intrinsic dispersion of the 
SN la peak magnitudes. Thus, for uncorrelated Gaussian errors, we expect X to follow 
a normal distribution with mean fix ~ and dispersion ax ~ 1- 

The angular correlation function of the magnitude residuals for a certain angular 
separation is defined as the expectation value of the product of the magnitude residuals 
of all SN pairs separated by that angle, 

C(0) = {X . X(0)) , (3) 
where the angular separation, 9, between two SNe (denoted 1 and 2) is given by 

cos(#) = sin(90°-5i) sin(90°-5 2 ) cos(ai-« 2 )+cos(90°-5i) cos(90°-5 2 ) -(4) 
Note that our definition of the correlation function is equivalent to the correlation 
coefficient for the original residuals, 5m, 

C(9) = ps m 8m(e) ■ (5) 
If C(9) > 0, the magnitude residuals are correlated, if C(9) =0 they are uncorrelated 
and if C(9) < they are anticorrelated. 

2.2. Detection limits for Gaussian magnitude residuals 

Under the assumption that the magnitude residuals follow a Gaussian distribution with 
a mean of zero and a variance a\ ~ 1, their product will follow a probability density 
distribution given by 

p(X • X) = ^K O^S\ , (6) 
™x V °x J 

where K is the modified Bessel function of the second kind and X ■ X denotes all 
possible products between the magnitude residuals. The variance of this probability 
density distribution is c\. x = °x- 

The correlation for a given angular separation is the expectation value of the 
product of the magnitude residuals for N p SN pairs, where N p is the number of unique 
pairs that can be defined from a set of iV SNe, 
N 2 - N 

= (?) 
The dispersion in this expectation value is given by 

*c = ^ = ^L. (8) 



N p VN P 



Since a x ~ 1, we can define the n-standard deviation (n-a) detection limit of the angular 
correlation to be 
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3. Analysing current data for angular correlations 

In this section, we will apply the methodology described above to two current SN la 
data sets in order to search for possible angular correlations in the magnitude residuals. 
We again stress the importance of using as homogeneous a data set as possible in order 
to minimize the impact of systematic effects connected to, e.g., the light-curve fitting 
method used when deriving the peak magnitudes. 

3.1. Data sets 

The first year Supernova Legacy Survey (SNLS) data set in Astier et al j3] (we will 
denote this data set as astier06) contains a total of 115 SNe la, of which 71 are high z 
SNLS SNe (z > 0.2) and 44 are nearby ones (z < 0.2). Figured] (left panel) shows the 
distribution of the SNe in the sky in galactic coordinates. Whereas the nearby SNe la 
(plus signs) are scattered across the entire sky, the SNLS SNe (triangles) are confined 
to four pencil-beam patches, each covering a one-square-degree area. These patches 
are referred to as D1-D4 and contain 14, 8, 31 and 18 SNe, respectively. The SALT 
light-curve model was used to fit all the SN la light curves in the data set. We add an 
intrinsic dispersion of <7 int = 0.13 to the uncertainties in the observed SN la magnitudes 
provided. 

The second data set is from Davis et al [31 HI El [34] (we denote this as davis07) 
and is a compilation of 192 SNe la consisting of 60 ESSENCE SNe, 57 SNLS SNe, 30 
HST 'gold' SNe and 45 nearby SNe. Figure [1] (right panel) shows the distribution of 
the SNe in the sky in galactic coordinates. Most of the low z SNe (plus signs) are 
the same as in astier06 and thus fairly equally distributed across the sky. The 147 
high z SNe (diamonds) are, on the other hand, very unevenly distributed. We use the 
peak magnitudes and uncertainties reported by davis07, which were obtained using the 
MCLS2k2 light-curve fitter on all photometric data. The uncertainties include both the 
observational and the intrinsic magnitude scatter. 

Since the two data sets have many SNe in common, e.g., most of the nearby SNe, 
the corresponding correlation functions will not be independent. Comparison of the 
results will provide a useful check on the impact of the light-curve fitter used in deriving 
the SN la peak magnitudes on the angular correlation function. 

3.2. Results 

For a given data set, we find the best fit cosmo logical parameters, and subtract 
the corresponding magnitudes, m fit , from the observed magnitudes, m b s , to get the 
residuals, 8m. From these, we subtract the mean residuayl, ns m , and divide by the 
uncertainty, a m , to obtain the normalized and scaled magnitude residuals, X. 

The values of the magnitude residuals and, consequently, the correlation function 
obviously depend to some degree on what cosmology we subtract. However, by varying 

| The value of [i$ m is very close to zero. 
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astier06 - 115 SNe davis07 - 192 SNe 




Figure 1. Distribution of SNe la in the sky in galactic coordinates. Left panel: 115 
SNe la in astier06, consisting of 71 SNLS SNe (triangles) and 44 low z SNe (plus signs). 
Note the pencil beam geometry of the four SNLS survey patches. Right panel: 192 
SNe la in davis07, with 147 high z SNe (triangles) and 45 low z SNe (plus signs). Many 
of the SNe are the same as in astier06. 

the cosmological parameters of the subtracted model, we have verified that the impact 
on the correlation function is very small. 

The correlation functions for the 115 SNe la in astier06 and 192 SNe la in davis07 
are presented in figure [2j The error bars are obtained using Monte Carlo simulations 
where we generate 2000 new data sets by adding random numbers from a Gaussian 
distribution with a mean of zero and a variance of a\ = 1 to the magnitude residuals 
of the SNe la. These data sets are then analysed to investigate the possible spread in 
the correlation function. The error bars represent the 68% confidence limit. The data 
is binned in such a way that each bin contains approximately the same number of SN 
pairs with an angular resolution of approximately 15°. Note that since each SN la is 
included in several bins, the data points will not be independent from each other, as 
seen clearly in the right panel of figure [21 where almost all the data are below the line. 
Thus, we do not expect about 1/3 of the points to lie more than la away from the zero 
line, as would be the case for uncorrelated data points. The horizontal bars indicate 
the range of each bin and the points are placed at the average angular separation in 
each bin. Both data sets are clearly consistent with the SN magnitude residuals being 
uncorrelated, but allowing for correlations at the ~ 10% and ~ 5% level for the astier06 
and davis07 samples respectively. 

3.3. High z 

We expect the effects from inhomogeneous DE to be more prominent at large distances 
and also systematic effects from, e.g., peculiar velocities to decrease with distance. 
Therefore, we check the high z SNe la in both astier06 and davis07 for correlations 
separately, i.e. we exclude the low z SNe la in the analysis. The results are presented 
in figure O The error bars represent the 68% confidence limit and are obtained using 
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Figure 2. Correlation functions for the 115 SNe la in astier06 and the 192 SNe la 
in davis07. The error bars represent the 68% confidence limit. Each bin contains 
approximately the same number of SN pairs. The horizontal bars indicate the range 
of each bin and the points are placed at the average angular separation in each bin. 
Note that, the data points are not independent, since each SN contributes to several 
bins. 



the same Monte Carlo simulations as above. 

The left panel in figure [3] shows the correlation function for the 71 SNLS SNe in 
astier06. Due to the pencil-beam geometry and the locations of the four patches of the 
SNLS survey, the SNe will be separated only by certain angles. We therefore do the 
binning at precisely these angles. The results are consistent with the SN magnitude 
residuals being uncorrelated. We note that there are differences in SN la brightness 
between the patches, e.g., there is a difference in the mean scaled magnitude residual 
between patch Dl and D4 of X pa 0.65 (where D4 is the brighter), corresponding to a 
difference of {5m) ~ 0.11 mag in the mean values of the original residuals. However, 
this is consistent with the expected statistical fluctuations at the 2a level. A search 
for correlations inside the patches, i.e. on sub-degree scales, yields a zero result for 
magnitude correlations also on small scales. 

The right panel in figure [3] shows the correlation function for the 147 high z SNe 
in davis07 adopting a resolution of 22.5°. Again, the results are consistent with the SN 
magnitude residuals being uncorrelated. 

3.4- Low z 

We also analyze the low z SNe, in the redshift range 0.015 < z < 0.125, separately. As 
noted, correlations between SNe la at low redshifts could be due to, e.g., correlations 
in their peculiar motions. However, since the SNe la in our sample are in the Hubble 
flow, we expect the contribution to the peak magnitute to be small. The results for the 
44 nearby SNe in astier06 are presented in figure H] (left panel). The data are binned in 
such a way that each bin contains approximately the same number of SN pairs with an 
angular resolution of approximately 15°. We detect an anticorrelation at the 2a level 
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Figure 3. Correlation functions for the high z SNe la in astier06 and davis07. The 
error bars represent the 68% confidence limit. Left panel: correlation function for the 
71 SNLS SNe in astier06. The bins are at the certain angular separations defined 
by the survey geometry. Right panel: correlation function for the 147 high z SNe in 
davis07. The binning is done uniformly. The horizontal bars indicate the range of each 
bin and the points are placed at the average angular separation in each bin. 



around 9 ~ 40°. The 45 low z SNe in davis07 (out of which 41 SNe are also in astier06) 
show the same anticorrelation. Thus, we can exclude that this is an artefact of the 
fitting technique used for determining the peak magnitudes of the SNe. However, other 
systematic uncertainties, which are not included in the error bars, could contribute to 
decrease the significance of this detection. 

In the right panel, we plot the scaled magnitude residuals, X, at every point in the 
sky by taking the weighted average of the residual for all SNe. The weight for each SNe 
is exponentially decreasing with respect to the distance to the point where the average 
is computed, using an exponential scale length of 15°. 

4. Future data sets 

As already noted, different SN surveys have very different survey geometries. The survey 
geometry affects the distribution of the angular separations which, in turn, affects on 
what scales and with what precision it is possible to detect possible correlations. We now 
consider the expected detection limit for SN la magnitude residual angular correlations 
for four different survey geometries, corresponding to the complete SDSS-II, SNLS and 
SNAP surveys as well as an all-sky survey. 

4.1. SDSS-II 

The SDSS-li§| (Sloan Digital Sky Survey) Supernova Survey's survey area, called 'Stripe 
82' (southern equatorial stripe), is a region 2.5 degrees wide along the celestial equator 



§ www.sdss.org 
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Figure 4. Left panel: correlation function for the 44 low z SNe in astier06. The error 
bars represent the 68% confidence limit. Each bin contains approximately the same 
number of SN pairs. The horizontal bars indicate the range of each bin and the points 
are placed at the average angular separation in each bin. Right panel: colour map 
of the scaled magnitude residuals for the 44 low z SNe in astier06 smoothed over an 
angular scale of 15°. 



(i.e. —1.25° < 5 < 1.25°), from roughly —60° < a < 60°. They aim of the survey is to 
find roughly 500 SNe la with 0.1 < z < 0.3 in this area out of which 300 SNe la are 
expected to be used to build the Hubble diagram. 

4.2. SNLS 

The Supernova Legacy SurveytJJ] (SNLS) is a pencil-beam survey, using four small square 
patches, each covering a one-square-degree area (see figured] in section [3Tl for the patch 
configuration). This is an example of a survey geometry that limits the detection of 
possible correlations to certain angular scales. One can either look for correlations 
inside the patches or between patches (six different combinations). Around 500 SNe la 
with redshifts up to z ~ 1 are expected to be used to build the final Hubble diagram. 

4.3. SNAP 

The survey geometry of the upcoming SNAlfl (SuperNova Acceleration Probe) satellite 
is not yet determined, but the plan is to cover a 7.5-square-degree field in both the 
northern and southern hemisphere. In total, SNAP is expected to find 2000 SNe la with 
z < 1.7. We assume that the patches have a square geometry. 

|| www.cfht.hawaii.edu/SNLS 
% snap . lbl . gov 
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44. All sky 

Pan-STARRsE] (Panoramic Survey Telescope and Rapid Response System) will as a by- 
product be able to detect a large number of SNe la across the entire sky. The proposed 
LSST0 (Large Synoptic Survey Telescope) will search for faint astronomical objects 
across the entire sky and detect large amounts of SNe la. Depending on the observational 
effort put into following up these detections in the future, it is not unreasonable to 
assume that this could result in up to 2000 SNe la being put on the Hubble diagram at 
low and intermediate z within the next decade. 



4-5. Results 

We generate random SN positions in the ranges defined by the survey geometries and 
compute the 1-a detection limits C\ im for each of the four surveys described above. 
In the results presented in figure El we have assumed that the magnitude residuals in 
each bin correlate in the same way, i.e., the error bars do not reflect the fact that the 
correlation function may vary over the size of the bin. 

The amplitudes and detailed shapes of the individual curves depend on the desired 
resolution. The lower the resolution, the larger the number of SN pairs in each bin and 
the lower the detection limit. We can scale the curves by making the crude estimate 
that the number of SN pairs in each bin scales with the bin size A8 as N p oc AO, such 
that Cii m decreases by a factor of \JK if the resolution is degraded by a factor of K. 

For the all-sky survey, there will be a peak at 90° in the distribution of the angular 
separations, with relatively few SN pairs at the smallest and largest separations. The 
all-sky survey thus probes correlations around = 90° most effectively. In the top left 
panel of figure El a resolution of 12° is assumed. The top curve is for a data set with 
1000 SNe and the bottom curve for 2000 SNe. 

For the SDSS-II, the angular separations are distributed with a maximum at small 
angles and then fall off for increasing separations. Beyond 120° it is not possible to 
probe possible correlations with this survey. In the top right panel of figure El the two 
curves are for a data set with 300 SNe, with an angular resolution in the top curve of 
10° and in the lower curve of 20°. The SDSS-II has its lowest detection limit at small 
angular scales. 

Due to the pencil-beam geometry of the SNLS, it is only possible to probe certain 
discrete angular separations in this survey. We generate 125 SNe in each of the patches 
D1-D4 described in section [3TTI for a total of 500 SNe. This means that there are equally 
many SN pairs in the bins, except for the first one, which probes correlations inside the 
patches. This bin contains more SN pairs and thus has the lowest detection limit as is 
evident from the lower left panel in figure El 

For the SNAP survey, we generate 1000 SNe in each patch, for a total of 2000 SNe. 
The distribution of the angular separations inside the patches resembles the shape for 

+ pan-starrs.ifa.hawaii.edu 
* www.lsst.org 
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Figure 5. 1-er detection limits for correlations in the magnitude residuals in different 
surveys. Top left panel: all-sky survey with 1000 SNe (top curve) and 2000 SNe 
(bottom curve). Top right panel: SDSS-II with 300 SNe for two different bin sizes. 
Bottom left panel: SNLS with 500 SNe. Note that the pencil-beam geometry only lets 
us probe certain discrete angular separations. Bottom right panel: SNAP with 2000 
SNe for 9 < 3° . The overplotted dashed line is the lensing correlation of equation (TTTj) 
with z s = 1. 



the all-sky survey but with a small tail at the high end. In figure [5] we have cut off the 
part of the distribution with 3° < 8 < Vl5° for readability reasons. 

4-6. Lensing correlations 

As noted in the introduction, physical effects such as peculiar motions and gravitational 
lensing will introduce correlations in the SN la data. Out of these, gravitational lensing 
will be the dominant contributor at z > 0.2 |35j. The covariance of the weak lensing 
convergence k for two sources at the same redshift z s can be estimated using [36] 

/ a \ -(n+2)/2 

(K.K(e))V^omaA 7 ^\Y°) ' (10) 
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where n is the spectral index of the power spectrum of density fluctuations. For small 
convergence (k <C 1) the magnitude residual correlation due to weak lensing is given by 

where we have assumed as = 0.8, Q m = 0.3 and n = — 1. It is evident that correlations 
from gravitational lensing will only be important at small angular scales where we 
expect DE inhomogeneities to be negligible. A tentative detection of the effects of 
gravitational lensing on SNe la in the GOODS fields was made in Jonsson et al [37] by 
cross-correlating the SN la residuals with the lensing properties of foreground galaxies. 
The dense sampling of the SNAP survey fields will allow us to detect the angular 
correlation of SN la residuals without using the information from the lensing galaxies. In 
the bottom right panel of figure the lensing correlation for a source redshift of z s — 1 
(approximately equal to the average SN redshift expected for SNAP) and a m = 0.2 is 
compared with the SNAP correlation detection. It is clear that SNAP will be able to see 
the lensing correlation with a resolution of < 0.3°; see also [38] for how this information 
can be used to probe cosmological parameters. We also note that with the increased 
statistics, it will be possible to probe the angular correlation function in separate redshift 
intervals, or to do a full three-dimensional correlation analysis. 



5. Toy model 

In this section, we assume a toy model angular correlation function, and illustrate how 
well we can detect such a correlation using future data expected for the SDSS-II and 
SNLS surveys described in section HI We let the magnitude residuals be the sum of an 
uncorrelated part, X ; , and a correlated part, X c , due to, e.g., inhomogeneous DE, 

X = Xi + X c . (12) 

Our toy correlation function has the form of a Gaussian, 

C(9)=Aexp{-[(9-e a )/9 c ] 2 }, (13) 

where the amplitude, A, is the variance of the correlated magnitude residuals, a\ c = A, 
8 a is the angle where the correlation has its maximum value, and 8 C is the angular scale 
of the variation in the correlation function. 

We generate random SN positions and uncorrelated intrinsic magnitude residuals 
with a mean of zero and a variance of a\. = 1 — A. The ratio of the dispersions of the 
intrinsic and the correlated magnitude residuals is 



^=V^- (14) 
ox V A 

Figure [6] shows the correlation function equation ffTB"]) with 9 C = 30° and 6* a = 0°. 
The error bars represent the 68% confidence limit obtained from equation (JH]) for 300 
SDSS-II SNe (triangles) and 500 SNLS SNe (diamonds). We have neglected the fact 
that the correlation function varies over the range of each bin. The binning is uniform 
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Figure 6. Our toy correlation function with 9 C = 30° and 9 a = 0° for 300 SDSS-II 
SNe (triangles) and 500 SNLS SNe (diamonds). The error bars represent the 68% 
confidence limit. Left panel: correlation function amplitude A = 0.1, corresponding 
to (JXi/vXc = 3. Right panel: A = 0.01, corresponding to <JxJ&x c ~ 10. The SDSS- 
II has uniform binning while the SNLS correlation is evaluated at the seven discrete 
angular separations defined by the survey geometry. 



for the SDSS-II. For the SNLS, we have used the discrete angular separations defined 
by the survey geometry. The left panel shows the correlation function with A = 0.1, 
corresponding to <Jx-Jvx c = 3. The correlation function is detected at high confidence 
using either data set. The right panel shows the correlation function with A = 0.01, 
corresponding to Ox-J<Jx c ~ 10. The SDSS-II data set is not sufficient for detecting the 
correlation in this case. For the SNLS, we expect to be able to make a weak detection 
of the correlation within the four survey patches. 

6. Conclusions 

The cosmological community is hard at work trying to constrain the behaviour of DE. 
The most pressing question is whether DE is a CC or something dynamical. A detection 
of temporal or spatial variations of DE would answer this question and refute the 
assertion of the CC being the dominant energy component in the universe. At the same 
time, there are alternative attempts to explain the apparent acceleration of the universe 
without invoking DE at all: instead as originating in the large scale inhomogeneities in 
the matter distribution. Irrespective of their origin, large scale inhomogeneities would 
manifest themselves as anisotropies in the observed magnitudes of SNe la. 

In this paper, we have devised a methodology for detecting angular correlations 
in SN la magnitude residuals. The methodology was applied to two recently found 
data sets (astier06 and davis07), neither of which show any signs of correlations at 
angular scales 0° < 9 < 180°. The uncertainties on the measured correlations, C(Q), 
are approximately 10% and 5% for the astier06 and davis07 data sets respectively, 
using an angular resolution of ~ 15°. Note, however, that the two data sets are not 
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independent, since they overlap partially for the high redshift sample, and include almost 
the same set of nearby SNe. The main difference between the data sets is that different 
assumptions on SN la properties, such as intrinsic colour and dust extinction in the 
host galaxy, are used when deriving the SN la peak magnitudes. Such assumptions are 
responsible for some of the most important systematic uncertainties in SN la cosmology. 
A comparison between the results for the two data sets allows us to keep the impact 
of these uncertainties on the measured correlation function under control. Because 
systematic uncertainties are typically different at low and high redshifts, and we expect 
effects from DE inhomogeneities to grow with redshift, we have applied the same analysis 
to the two sub-samples separately. At low z, we found an anticorrelation at the 2a level 
at angular scales 9 ~ 40°. Due to the weak statistical significance of this detection, and 
the low redshifts, we cannot draw any conclusions about the possible implications for 
DE inhomogeneity. It is worth noticing that the nearby SNe analysed are a collection of 
very inhomogeneous data samples, collected by different observers at various telescopes 
and thus, there could be systematic uncertainties involved which have not been included 
in the error bars (see, e.g., discussion in [39] ). We note that we see no signs of anisotropy 
on angular scales 9 ~ 180°, corresponding to the north/south asymmetry detected in 
|33j. Our data sample, however, is only partially overlapping with that work, since 
we have only included SNe in the Hubble flow (z > 0.015) in our analysis in order to 
mitigate the effects from correlations arising from peculiar motions. 

Gravitational lensing correlations will only be significant on small scales, where we 
expect the effects from DE clustering to be negligible. Using data from the proposed 
SNAP satellite, we should be able to detect the lensing correlations in SN la magnitudes 
with an angular resolution of < 0.3°. 

Correlations induced by intervening galactic dust will also mainly take place on very 
small angular scales, whereas correlations from dust in the Milky Way could potentially 
be a problem for an all-sky survey. Most surveys, however, observe in directions in the 
sky where the Milky Way dust extinction is well measured and understood. 

More SN data are needed if we are to detect presumably weak correlations and 
better constrain inhomogeneous/anisotropic models. Our simulations of future data 
illustrate how the survey geometries govern on what scales, and to what precision, we 
will be able to detect possible angular correlations. We find that using data from the 
soon to be completed SDSS-II and SNLS surveys, we will typically be able to detect 
correlations in the magnitude residuals at the per cent level. 

Any claim of a deviation of DE properties from those of a CC will be subject 
to intense scrutiny from the cosmological community, and will need to be backed up 
by independent evidence in order to be generally accepted. It is therefore of utmost 
importance to pursue the study of DE along multiple paths. Although observationally 
demanding, the search for spatial variations in DE properties is a useful complement 
to studies aiming at constraining the time evolution of DE. Specifically, for SNe la, 
we expect systematic effects connected to the time evolution of DE, such as temporal 
evolution of SN la properties, to be of less importance when studying spatial variations 



Probing dark energy inhomogeneities with supernovae 



15 



of DE. Since systematic effects constitute the limiting factor for future SN la surveys, 
we expect interesting results for the spatial correlations of SN magnitudes in the near 
future, irrespective of whether such correlations are found or not. 
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